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We identify a two-parameter family of excited states within slow-roll inflation for which either 
the corrections to the two-point function or the characteristic signatures of excited states in the 
three-point function - i.e. the enhancement for the flattened momenta conflgurations- are absent. 
These excited states may nonetheless violate the adiabaticity condition maximally. We dub these 
initial states of inflation calm excited states. We show that these two sets do not intersect, i.e., 
those that leave the power-spectrum invariant can be distinguished from their bispectra, and vice 
versa. The same set of calm excited states that leave the two-point function invariant for slow-roll 
inflation, do the same task for DBl inflation. However, at the level of three-point function, the calm 
excited states whose flattened conflguration signature is absent for slow-roll inflation, will lead to an 
enhancement for DBI inflation generally, although the signature is smaller than what suggested by 
earlier analysis. This example also illustrates that imposing the Wronskian condition is important 
for obtaining a correct estimate of the non-Gaussian signatures. 
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I. INTRODUCTION 

Primordial non-Gaussianity is among the most promising probe of our early universe [l[. Although current cos- 
mological data is consistent with the primordial density fluctuations being Gaussian, any observed departure from a 
Gaussian spectrum will likely point us to some interesting microphysics. In the simplest vanilla version of inflation, 
namely, slow-roll inflation with a Bunch-Davies intial state, the primordial non-Gaussianity produced was shown to 
be suppressed by the slow-roll parameters 0,11] and is thus unobservably small. Nonetheless, well motivated deviation 
from this standard picture can lead to a detectable signal in current and upcoming experiments. For example, even 
within the context of single-field inflation, non-canonical kinetic terms can generate large bispectra (3-point function) 
of the equilateral shap es 11 , [Bl ; a non Bunch-Davies initial state of inflation can boost the folded/flattened momentum 
configuration 0, 0, Il0l|: and features in the Lagrangian can result in a large bispectrum with oscillatory running 
A combination of the above effects realized as resonant folded non-Gaussianity was considered in [1^. Mul- 
tifield inflationary models open up even further possibilities, such as curvatons f)3\, turning [l"5l - l22j or bifurcating 
[Ullll trajectories, thermal effects [Ull^l and etc [H]. 

Among the various possibilities, the choice of initial state of inflation is perhaps the least understood from a 
microphysical point of view. While the other scenarios mentioned above can be analyzed within well deflned field 
theoretical models, the choice of initial conditions for inflation require a detailed understanding of short distance 
physics. The fact that the universe may have inflated much more than needed to solve the flatness and horizon 
problems of standard cosmology implies that the CMB scales today may have been not far from the Planck scale 
at the beginning of infiation. This is sometimes known as the "trans-Planckian problem" of inflationary cosmology 
[28j . though we hasten to stress here that this problem is not specific to Planck scale physics. We will hereforth 
refer to the Trans-Planckian problem as the question of how to formulate infiation when there exist a short distance 
scale M (e.g., the string scale Mg) not far from the Hubble scale of inflation. Our understanding of physics at such 
high energies is rather limited at present. One may however view this UV problem of infiation as an opportunity, 
as CMB and other cosmological measurements can provide us a way to probe high scale physics. The short distance 
modifications of inflation can affect the dynamics of the primordial perturbations [l^, [s^l as effective field theory 
arguments would suggest [3l|, [s^ , or their initial conditions [s^l • While a UV completion of infiation that allows us 
to unambiguously determine its initial state is not yet available, one often resorts to a phenomenological approach. 
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In [33|, a class of initial states was motivated based on the notion of vacuum (i.e., minimizing particle production) 
in a theory with a UV cutoff. The computed power spectrum turns out to be different from that obtained with the 
Bunch-Davies vacuum. Subsequent work in formulating such initial state effects in terms of boundary effective field 
theory [sJ] further exemplified this fact. Thus, it is generally expected to see modification of the power spectrum 
with any non Bunch-Davies initial state though the deviation may not always be large enough to be observed. 

Having an excited state as the initial condition for inflation can be justified as a pre-infiationary phase could have 
excited the infiaton fields from their vacuum (35l.l36j or some non-adiabatic process, at the time when the energy scale 
of inflation was around Af , could excite the vacuum to an excited state 22, 37]. The effect of such excited states 
are usually discussed in terms of modulated oscillations on the power spectrum. However, as originally pointed out 
in Q, the effect of such excited states may be more conspicuous at the level of three-point function probes, especially 
for slow-roll inflationary models with higher derivative correction terms [1, 0, Q or in the context of DBI inflation 

In this note, we address a related but different question: can such initial state effects be hidden? We identify excited 
states that can deviate from the Bunch-Davies vacuum as much as it is allowed by the constraints from backreaction, 
but do not leave any signature at the level of the two point function. We call these excited states "calm excited 
states". We show that all these excited states are revealed at the level of the three-point function through their 
enhanced contribution in the flattened momentum configurations. Our analysis also suggests that an approximation 
which was made in previous studies [1, , in which the first Bogolyubov coefficient was set to its Bunch-Davies 

value, may lead to an inaccurate estimate of such enhancement. We argue the same set of calm excited states would 
leave the power-spectrum invariant in DBI infiation. For slow-roll inflation, one can also identify excited states that 
again maximally violate the adiabaticity condition but the enhancement for flattened configurations is identically zero. 
However, in DBI inflation, it it is not possible for an excited state to give vanishing enhancement to the flattened 
conflguration. 

The outline of this short note is as follows: in the first part we construct initial excited states that leave the two-point 
function of slow-roll inflation intact. We then compute the three-point function for a general excited state without 
resorting to any approximation for the Bogolyubov coefficients by imposing properly the Wronskian condition. We 
further identify the calm excited states in DBI infiation. Finally we show that in DBI infiation any excited state 
that maximally violates the adiabaticity condition could be observed by non-gaussianity probes unless the scale of 
new physics, M, is higher than 10~^ Mpi. Explicit expressions for the bispectrum with general excited state in DBI 
inflation are relegated to the appendix. We end with a brief conclusion. 



II. SLOW-ROLL INFLATION 



A. Calm Excited States with no Modification to the Power Spectrum 



1. Power Spectrum 

We would like to know if one can modify the initial conditions for perturbations that correspond to excited states, 
but still leave the power spectrum intact. Let us assume for simplicity that the inflationary background is de-Sitter 
space, a{r]) = The perturbations in such a background satisfy the following differential equation 

< + ik' - ^)uk = 0. (1) 

The full solution to this differential equation is 

u.irj) = (-ry)i/2i/« (-M + (-,7) ^^//(J) (_fc^). (2) 

The term proportional to Hy]^{—kr]) behaves like positive-frequency wave in the inflnite past and the other term 
behaves like the negative frequency one. The Wronskian condition, u*u' — uu'* = —i, will result in the following 
relation between ak and /3fe 



\ak\'-Wk\' = l. (3) 



The power spectrum for the general solution ([2]) is 



Ps = tei^ (4) 
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We would like to know if it is possible to obtain an unmodified scale-invariant power spectrum, Ps = starting 
with non-zero Uk and i.e. if it is possible that 



W-Pk?^\. (5) 



afc and Pk are in general complex variables. We write them as 



au ^ xi+ ix2, (6) 
/3fc = i/i + iy2- (7) 



Eqs. ([5]) and ^ could then be written as 

^i+Vi- 2a;iyi + + - 2x2y2 



4' 

xl+xj-yj-y^ = J. (8) 



From the above two equations one can solve for xi 



^^^yJ±A.^, (9) 

Vi Vi 

Plugging this back to the second equation of ([8]) and solving for X2 , one obtains a second order equation for X2 whose 
solutions are: 



xf=y2- '^^'^ . (11) 



One can then solve for xi using pO)) and (fTTj) and respectively and find two solutions for Xi: 



yi - y2 sign(?/i )—===, (12) 
2\/yJ + vi 



a;f ^ = 2/1 + y2 sign(?/i)— ^^==. (13) 
2Vyi + 2/2 

Let us mention here some illustrative examples of calm excited states: Pk can be chosen such that the number of 
particle that are produced during inflation do not backreact and spoil inflation. In principle, the real and imaginary 
parts of j3k do not have to be related, but for simplicity, let's consider a first example where they do: 

/3fe=a(l + z). (14) 

From the bounds on backreaction, a is expected to be smaller than Q 



a<mm{Ve^^,V^^^}. (15) 

M is the momentum cut-off, below which the effective field theory description is valid. The bound on /3k comes from 
the fact that the energy density of the nearly massless quanta of the inflaton has to be less than MpH^. Also the 
produced energy should not violate the slow-roll condition, hence the appearance of the first and second slow-roll 
parameters, e and rj' . These two considerations will lead to the above bound on a. In obtaining the results, one has 
to assume that /3fc — >■ for k > M(a(ryo))i where 770 is the conformal time the physical momentum of the comoving 
mode becomes equal to the cutoff of the effective field theory, i.e. k/a{rjQ) ~ AI [6] 38]. From the above consideration, 
one can see that a is not necessarily very small. For example, for an infiationary potential with V = ^to^0^, with 
H ~ 3.6 X 10~^ and e ~ 77' ~ 0.01, for M = 107J, a ~ 1.36 in the above example. Note that earlier analysis in 
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the literature which constraint M to be greater than lOOiJ 0, [1] do not apply to the the above example as the calm 
excited states are special in that they do not leave any signature in the power spectrum. 

For the above choice of /3fc, eg. (fT4| . the first Bogolyubov coefficient could be obtained from eqs. (fTOl) - (fT3l) to have 
a calm excited state: 

and 

Starting with either pair of {ak^^\l3k) and (afc^^\/3fc), one obtains the same power spectrum as that for the Bunch- 
Davies initial state. 

Another example which maximally deviates from the Bunch-Davies vacuum but still satisfies the bound from the 
backreaction is given by: 

h = CT. (18) 

The Wronskian condition then implies the following two solutions for a^: 

afeW=a + ^^, (19) 

and 

a,<•'^=a^^^. (20) 

The fact that for this case the imaginary part of (3k is zero will be of interest in the next section. As we will see, 
following |6| , one would obtain zero enhancement for flattened configurations for the above excited state, whereas the 
actual enhancement is nonzero if the ak coefficient satisfies the Wronskian condition ([3]) properly. 



2. Three-point function 

Let us calculate the three-point function for the above excited states to see if they modify the three-point function, 
and if so, how large the effect will be. One can calculate the Wightman function for the solution ^ 

G.>(.,.')-5^4^ (21) 
02 a{r]) a{ri') 

The three-point function could be derived from the Wightman function through the following integral 0: 

+ permutations + c.c. (22) 

To compute the three-point function, the first argument of the Wightman function has to be set to zero, which 
corresponds to the moment the mode is outside the horizon, and then differentiate it with respect to the second 
argument. We obtain 

a{ri)dnG^ (0, 77)) = ^ (23) 

Plugging this result back to equation (|22|) . one obtains: 

3 /-r \ J=3 



(C-.C^Cr>) = (T ^ exp(zfc.,) - ^^^^^^ - S^^:' V exp(-z^,) 

-\- c.c. + permutations (24) 



\ 

_3exp(-iA:t?7) 
+ Pfc 
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where ijo is the moment at which the physical momentum becomes equal to the physical cutoff, i.e. r]o 
the sum of the magnitudes of the momenta kt — ki + k2 + k^ and kj ~ kt — 2k j. Upon integration, one obtains; 



M 
Hk ■ 



kt is 



where 



(1 - cos(fct7;o))^ , .sinktrio 

; Ol + I ; G2 



(1 - cos(fcj-)?7o) 



i=i 



sin kj r/o 



(— J- 



fc| 

(25) 



C4 = [ak - f^kfim^'^k -amf) + [c^ -'^fiakl3k^ ~ ak'^Pk) 



(26) 
(27) 
(28) 
(29) 



As it was noticed originally in the fact that kj in the denominator of d^S]) vanishes, enhances the contribution of 
these flattened configurations with respect to the Bunch-Davies vacuum. The divergence in the kj ^ is removed by 



considering that the cut-off is only relevant for kj 
one obtains: 



Thus for the configurations in which one of the kj vanishes, 



ka' 



h'^kik2ks 



(1 - cos(fct??o . smktrjQ 

; 61 + I ; 62 + iCiTjo 

Kt kt 



( i,2 + l2 + l2 ) 
^1 ^2 '^3 



(30) 



The first two terms will change the normalization and shape dependence of the three-point function. The last term is 
the enhancement which is important for the flattened configurations. The relative enhancement factor with respect 
to the standard result is: 



iCiktrjo — iCi 



k.i=0 



a{rio)H 



(31) 



The result is enhanced by the ratio of the physical momentum to the Hubble scale at the beginning of inflation. The 
factor C4 is what has been already approximated with 5(/3fc) in Q, by setting Uk = 46]. Here, we do not resort 
to this approximation. C4 could be factorized as 



C4 = i\akf - |/3fe|')(|afe - Pkf){akl3k-c^Pk) 
Using the Wronskian condition ([3]) and the condition ([5)), we can rewrite it as 

C4 = — ^(a^^fc) 
o 

- 2 

ITT 



(32) 

(33) 
(34) 



Note that setting ak — we recover that C4 is proportional to 3(/3fc), as stated in [(>]. Using the pair of expressions 
((ITni),([ni)) or (dni), dlSl)), one respectively flnds: 




yilyilv^ , Vl sign(yi)0F 



vl 



^'-^-^Vi+vl sign(2/i) 



(35) 



From the last expression for C4 one can see that in contrast to [6| , the enhancement for the flattened configuration is 
never zero for the calm excited states, even when the second Bogolyubov coefficient is completely real. Therefore the 
flattened configurations are enhanced even in the case of (fT8|) . The enhancement is proportional to 
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(36) 



This illustrates the importance of the Wronskian condition in getting the correct estimation of non-Gaussianities. 
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B. Excited States with No Enhancement for the Flattened Configurations in Slow-roll Inflation 



In the previous sections, we investigated whether there exist calm excited states that leave the two-point function 
unmodified and noticed that all these excited states will lead to an enhancement of the flattened configuration. Now 
that we are equipped with the general expression for the enhancement factor, let us see if we can also identify a 
family of excited states that do not lead to any enhancement of the bi-spectrum for the flattened configurations. Such 
excited states still satisfy the Wronskian condition, eq. ([3]). However instead of equation ([5|), which was necessary for 
obtaining an intact two-point function, we would like to assume that the enhancement for the flattened configurations 
is absent, namely we have 

C4 = (lafep - \/3k\^){\ak - l3k\^){aA-c^(3k) = (37) 

which can be simplified further using the Wronskian condition and writing the expressions in terms of real and 
imaginary parts of ak and 

C4 = ^(l"fe - Pk\'^){ak/3k - afc/3fc) = -■i^{y2Xi - X2yi){yl + - 2x2y2 + xj - 2yixi + yl). (38) 

Finding such nontrivial excited states is therefore contingent upon whether non-trivial solutions for the following pair 
of equations exist: 

x\+xl-yl~yl = J (39) 

{y2Xi - X2yi){yl + xl - 2x2y2 + xl - 2yixi + yl) = (40) 
The equation (liU]). by itself constitutes of two equations: 

y2Xi - X2yi = (41) 

yj+xj- 2x22/2 +xl- 2yixi + yl = (42) 

Solving p9l) and (j42p . one would obtain solutions for xi and X2 (in terms of j/i and 2/2) that are not real. Therefore, 
these pairs of equations would not lead to any meaningful results. Solving the other pair of equations p9p and (j4ip 
together yield the following results for xi and X2 

X2 = y2[ 2,2 (43) 



j + yl + yi 



These two equations determine a two-parameter family of solutions for which the enhancement from flattened config- 
urations are absent. Two example of such states are 



Pk 
ak 

and 



- CT(l + i) (45) 

- (^+aY/^+z(^ + aY/^ (46) 



> 1/2 



(47) 

«/c = + a') (48) 

One can work out the power spectrum for this class of calm excited states and see that the power spectrum takes 
the form 



l 7T~^^7T + A\(3k\'\l3k\+8\Pk\' \ 
^ ~ (2^)2e I - I (^^^ 

For very small values of \lik\, the correction to the power spectrum varies linearly with \ j3k\ for such excited states 

Ps ^^(1-^IA-l) (50) 
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III. DBI INFLATION 

The perturbations in the DBI inflation [s^, satisfy the following equation of motion: 

v'l + clk'vk - —vk = (51) 
z 

where Vk is related to the Fourier transform of the gauge- invariant perturbations of the comoving hypersurface, Ufc, 
through the following relation 



Vk = zCk, z = (52) 

Above, Cs is the speed of sound and e = —H/H^ is the first generalized slow-roll parameter. In de-Sitter space, 
a = — Tj^, the solution to differential equation (|5T|) could be expressed in terms of Hankel functions of the first and 
second kind: 

Vkiv) = akV^jH^yli-kcsTi) + l3kA/~^Hf'^^{-kcsri) (53) 

The Wronskian condition for the general solution of (|53p reduces to the previously derived equation, eq.Q. The 
two-point function for the general solution (j53|) is 

"^-"-^ 

Requiring the primordial power spectrum be unmodified, one again obtains equation ([S]). Since the two equations 
that determine the calm excited states (that leave the two-point function invariant) for DBI inflation are the same as 
that for slow-toll inflation, one can conclude that this class of calm excited states are common among both slow-roll 
and DBI inflationary scenarios. 

It is also useful to find the enhancement of flattened configurations for such calm excited states in DBI inflation. 
The three-point function in the limit that kj — > is calculated in the appendix. The leading terms in the limit where 
one has an flattened configuration is: 

A(C(i, ki)C(i, k2)C(t, k3))|fc^^o = £2f3{h){ktvof + £3f2{h){ktmf + £2fi{h)ktr^o, (55) 

with the explicit expressions for £2 and £3 given in the appendix. The factors fj{ki), j = 1 • • • 3, could also be obtained 
from the results in the appendix. For the calm excited states that leave the two-point function intact, one has: 

3 

£2 = T^sign(yi)|/3fe| (56) 

£3 = -l^f^kl'. (57) 

From the above dependence of £i on one sees that even assuming that the excited states maximally violate the 
adiabaticity condition, i.e. when condition (|15p is saturated, the quadratic and linear order (in fctryo) corrections will 
be negligible. This is because any probe of non-gaussianity through the CMB will be two-dimensional and such a 
projection from three dimensions to two, will cause one to lose a factor of \ktriQ\ in the enhancement Thus the 
real enhancement in this case comes from the cubic term. The enhancement is always nonzero for these states, as £2 
is proportional to \/3k\- 



A. Absence of Excited States with No Enhancement for the Flattened Configurations in DBI Inflation 

One may wonder if in DBI infiation it is possible to identify excited states for which the enhancement of flattened 
configurations is identically zero. It is easy to verify that there are no such excited states in DBI inflation, since 
the factors £2 and £3 could not be set both to zero simultaneously. One can identify at most the excited states 
for which the leading contribution to the enhancement in the non-gaussianity is zero, i.e. £2 = 0. Then the term 
linearly proportional to fctTyo will also be automatically zero. Nonetheless, the quadratic term proportional to £3 will 
survive. Since £2 oc C4 = 0, the excited state in DBI inflation for which the leading enhancement for the flattened 
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configurations is zero, are the same excited states in slow-roll inflation for which there is no enhanced nongaussianity 
for the flattened configurations. The coefficient of the quadratic term for such states is 



£3 = + 2|/3fe|V4|/3fcP+7r-^|/3feP + (58) 

For small values of \(3k\, £z could be expanded as: 



Thus after taking into account that a two-dimensional projection of the three-point function will suppress the en- 
hancement by a factor of fcjTyoj the observable enhancement is 



A(C(t,ki)C(i,k2)C(t,k3)) 



(C(i,ki)C(t,k2)C(i,k3)) 



^ \lik\km (60) 

>0 



The factor fct?7o is the ratio of physical momentum at the beginning of inflation to the Hubble scale which may be as 
large as M/ H. Assuming that |/3/c| takes its maximum value allowed by the constraint in (|15p. the enhancement will 
be proportional to 



A(C(t,ki)C(t,k2)C(t,k3)) 



^ - (61) 



(C(i,ki)C(t,k2)C(t,k3)) 

Thus, the observable enhancement is only negligible if 

M > min{ V^}Mp (62) 

For an inflationary scenario with e c^i rj' ^ 10~^, the enhancement from the flattened configuration is absent if 

A/>10"2^#P, (63) 

Otherwise, the enhancement from the flattened configuration can be considerable. 

Here, we pause to compare the enhancement we found with the results of some previous works. In fio'l, the coefficient 
Qffc was simply set to Uk — throughout the calculation \^f\ without regards to the Wronskian condition. We believe 
this leads to an overestimate of the the enhancement factor for the flattened configurations to be proportional to (fcjTyo)'^ 
(before the two-dimensional projection). Our analysis shows that, not only but ak is important in estimating the 
effect of excited states in the three-point function. 



IV. CONCLUSION 



It is often expected that any deviation from the Bunch-Davies initial state of inflation would lead to a modification of 
the power spectrum. Here we showed that this expectation is not necessary true: we identified a two-parameter family 
of initial states in slow-roll and DBI inflation for which there is no modification to the two-point function. If it were not 
the backreaction constraint, the deviation could be arbitrarily large. All such excited states are distinguishable from 
their three-point functions as they lead to an enhanced non-Gaussianity in the fiattened configurations. A correct 
estimation of such enhancement requires one to take into account the Wronskian condition which has so far been 
neglected in some earlier analysis. Taking into account this condition, we identified a class of excited states in slow- 
roll infiation that maximally violate the adiabaticity condition for which the enhancement for flattened configurations 
is identically zero. In DBI inflation all excited states would lead to an enhancement for flattened configurations. 
This shows that for slow-roll infiation, a combination of two and three-point function measurements are necessary to 
detect the initial state. However in DBI inflation, the three-point function by itself can reveal the nature of the initial 
states for inflation. It would be interesting to see if such calm excited states could be motivated from a microphysical 
point of view, and whether additional features may arise in concrete models/scenarios. To present explicit expressions 
in this paper, we considered the simple case in which the the Bogolyubov coefficients are scale-independent. More 
realistic situations may involve coefficients in which they are scale-dependent but the results will be highly model 
dependent. Having a concrete model in hand, however, one can easily generalize our computations to incorporate 
the scale dependence of the Bogolyubov coefficients. Also, it would be of great interest to extend our current study 
to higher point functions such as the trispectra especially for models with non-canonical kinetic terms where the 
non-Gaussianities generated can be observably large ,41.-i44] . 



Appendix: Three-point Function for General Excited States in DBI Inflation 

For a general single field inflationary scenario with the following Lagrangian: 

S=IJ d^^V^ [MliR + 2P{X, cP)] , 
the interaction Hamiltonian for the curvature perturbations of the comoving hypersurface, C, to O(e^) is [1 
H;nt{t) - - /dM-«'(S(l-^)+2A)^ + 4^(e-3 + 3aCC' 

J Cg £1 Cg 

+ ^(e - 2s + 1 - cDmf ~ 2a4c(aC)(5x)} • 
(j) above is the inflaton field and X = —}^g^^'^dfj,(j)du4'- X satisfies the following equation: 



2,, 2 



ci 



Also 



X + r XX = — 5~ 



A = X'^P.xx + %X^P.xxx ■ 

It is useful to define the "speed of sound" Cg as 

2_dP _ P,x 



dE Px + 2XP^xx 
and some "slow variation parameters" as in standard slow roll infiation 

H XP,x 



The expectation value of the three-point function in the interaction picture is given by 

(C(t,ki)C(i,k2)C(t,k3)) rdi'([C(t,ki)C(t,k2)C(t,k3),iI,„t(t')]) , 

Below, we will find the contribution of each term in the interaction Hamiltonian: 

Contribution from C,^ term could be calculated in the following way. We denote kt = fci + ^2 + fcs, and find 

Ti = -z(c^-i + ^)^^o,kiMo,k2Mo,k3) 



d^^^^^") ^ (}^k.)+c.c. 



Using the following equations 

- I3k) 



u 



(0,k) = - 



du*{rj,'ki) iH'i]y^kcl{—ak exp{ikcsrj)) + (5k exp{—ikcsri)) 



drj y/ne 
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Ti will become: 
96 

Ti 



3 



5ft I i{ak - Pk) {oik I V exp{iktCsri) - jSk r] exp(-i/ctCsry)) 



+ y^Sft N(afc - /3fc)^(-afc^/3fc / rj"^ exp{ikjCsri) + (3k a]^ 

j—l V Jrio Jrio 



rf- cxp{-ikjCsri)) 



+ perm. 



(76) 



Upon integration yields: 



cfe2A;ifc2fc3 



(-24fi + cos(6)(24fi - 72Ctf2 - 12^1^') + sin(et)(72f2 + 24^1^* - 36f2^t )) 



3 I 

+ E Ti (24^3 - cos(^j)(-24£:3 + 24£:2e,- + 12£:3^|) + sin(^,)(-24£:2 - 24£:2e,- + 12£:3^|)) 



i=l 

where = ktCsTjo, = kjCsfjo and £i, i = 1..3 are as follows: 

£i = 2ak^af' - Sak'a^'^Pk + 3afc^/3fe^ - iak^OkPk + '^cukOk^k^k - ^otkPk^ h + ^oik'^Pk - ^oikhPk + '^Pk^ l^k 
£i = i{ak- lik)ip^- fik){-c^Pk + akfik) 

^3 = oikOk'^Pk + ak^f3k'^ + ak^OkPk - GakOkPkPk +akPk^Pk + ak^Pk + ockPkPk 
In the limit where one of the kj 's approaches zero, Ti would behave like: 

f - lim T - ^'^^'^^^ (f - 1 + Cs) 3 3 

For the contribution of term, one finds 

T2 = i4(e-3 + 3c^)u(0,ki)w(0,k2)u(0,k3) / a^dr] 

X 2K(,, k,) '^-^ '^^^ + sym)(2.)3^3(^ ^ 

Plugging the expressions for u(0, k) and j^'*^^ into the integral yields 
32if4(e-3 + 3c2 



(77) 



(78) 



(79) 



(80) 



-3? 



{ak - Pkf / (a^^(A;ic,?7 + i)e*'='"='' + pk {kic.r] - i)e 



(^e'^kfk2k3 



-iktCsT] 



^ i^-af PkikiCsT] + i)t 



akPk {kiCsTj - i)e 



-ikj CsT] 



j=2 



+ perm. 



(81) 



Performing the integration, one finds: 

^ ^ 7r^f/^i£_3 + 3c|) ^ 
cie^kik2k3 



+ (4.3^ + cos(6) f-4^3%# + ^1 + sin(a) (-^£2%- ''^^^ 



1 1 

HI2 



(Ml 



1 1 

+ E U£: 

j=2 \ 



;^+cos(^,)^4f3^pp-^^ 



+ sin(^,) 4£; 



perm. 



(82) 
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where 



1 ^ 1 

- = Y- 



i=l 

3 



- = T- 

7? - 1,2 



In the limit kj 0, we have: 

= Um Ta 



27r(-3 + 3c2 + e) 



1=1 * 



-^c^ryo + ^3 ( — 



2 2 



cfe2fciA;2A;3 
For the C(^C)^ term, one obtains: 

2ie /■° 
^3 = — (e-2s + l-c^)M(0,ki)u(0,k2)M(0,k3) / a^dr] 
c-s J no 

X ((k2.k3)w*(j?, ki)«*(r?, k2)w*(j?, ks) + sym)(27r)353(^ + c.c. 

i 

Plugging the expressions for the relevant parameters into the above, one obtains: 



(83) 

(84) 
(85) 

(86) 



(87) 



32g4(e-2s + l-c^)(k2.k3) 



5R 



{ak 



af{kiCsri + i){k2CsTi + i){kiCsr] + i)e*'=*'=»'' 



+ Jk [kiCsT] - i){k2Csr] - i){ksCsV - i)e~'''"''''' 

+ (ok^PkikiCsV + i){kiCs'n + iKhCsV - i)e'''^''''' + a^fikic^r] - i){k2CsV - i){k3CsV + i)e'*^^'=^'' + perm.))) 
+perm. (88) 
Upon integration, one finds: 



73 



t:H^{-1 + cl-e + 2s) / 2(fc^ + fc| + /c|) 



-3^2 cos(^( ) + £i sin(^( ) [-£2 cos(^j ) + £3 sm{^j ) 



+ perm. 



{kl + kl + kl) f 1 1 



(fcifcafcs)^ \ktli kl 

2 2 2 ^ 
~|~2^/t-|^ ~|~ ~\~ 



— + -2 (-3^2 sin(6) +£i{l- cos(^t)) 



{kik2k3fki 



(£:isin(^t) + 3^:2(1 -cos(6))) 6 



2 , ,2 , , 2^ I kliki + k2) + ksikj + kl) - kik2{ki + k2) 



{krk2k3fkl 



(£3(1 - cos(^3)) - £2 sin(C3)) + perm, j 



+2{kl + kl + kl) 



^ {k,k2k3Ykl ^-'^ + ^^"^^^^ V ^ . 



(89) 
(90) 



The contribution of the term in the first line is divergent when the upper limit of integral goes to zero. To make 
sense of it one should consider that the actual observation is made at the time of last scattering which interchanges 
the upper limit of integration with the corresponding conformal time at the last scattering surface, rf . In the limit 



of kj —^0, we have: 



T3 = lim T3 



27rif4(-l + c2-e + 2s)(k2.k3) 



{kik2k3)^cie^ 



[-£3^1^2^:30^770 + 252(^3 + ^2 - k2k3)cs'no\ + perm. (91) 



Finally the contribution of C(^C)(^X) will lead to the following integral 



Ta = - 



2ie-' 



u(0,ki)7i(0,k2)u(0,k3) / a^dr] 



X (2(k2.k3) 



(iM*(7?,ki) 



M*(?7,k2) 



du*{ri,'k3) 
dr] 



+ perm.)(27r)^^^(^ kj) + c.c. 



(92) 
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Plugging the expressions for the relevant parameters into the above integral, one obtains: 



Ta 



- {-^Wkik2CsV + i)e'^'''^''-am^ik2CsV-i)e'~'^'^'' 

J = {h3} 

-perm. 



Integrating yields 

87rij4(k2.k3) / 1 



£i{k2+kt) 



£i{k2 + kt) 



y ^ 



J = {1.3} 



£a{k2 + kj) 



£^{k2 + k 



k2 



a^aCt +sin(6) 



?>£2{k2 + kt) 



- £2ij + sm{^j 



k2 

£2(^2 + kj) 



^36 



1 



£a{k2 - h) 
k2 



+ cos(6) 



£3ik2 



k2 



£2^2 +sm{^2) 



^2(^2 - h) 



£3^2 



+perm. 

In the limit where fcj — >■ 0, we have: 

kj.kj+2 , kj-kj+i kj+i.kj+2 



r4 = Jmi r4 = 4 ■ ■ , 



(93) 



(94) 



(95) 



^3 



:jKj+2 



kj+ikj+2 



kj.kj+2 ^ kj-kj+i kj+i.kj+2 



'«j+2 



^2 1-2 



(96) 
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